Synchrotron radiation in a chromo-magnetic field 
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We study the generalization of QED synchrotron radiation to the QCD case with a chromomag- 
netic field using the Schwinger et al source method. It is shown that the QED case can be obtained 
as a special limit. The comparison with the path integral approach of Zakharov has shown consistent 
results. 
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I. INTRODUCTION 

It is believed that during and after relativistic heavy 
ion collisions strong chromomagnetic field will form that 
can be treated as a classical background. Numerical so- 
lutions P, [|| indicate that "just" after collision trans- 
verse color electric and color magnetic fields change sud- 
denly from being transverse in the initial state, in the 
so-called Color-Glass-Condensate state, to being longitu- 
dinal. The latter are called glasma flux tubes. Transverse 
fields then rise, and at some stage after the collision the 
transverse and longitudinal components of color electric 
and color magnetic fields reach a "steady" comparable 
values. In this rich environment a fast parton escaping 
from the collision would feel the effect of such a field. 
Synchrotron and Cerenkov effects are important physi- 
cal phenomena to be studied, either as an energy loss 
mechanism or as a coherent gluon radiation process. In 
this context the Cerenkov radiation was considered in 
Q . Syncrotron radiation was analyzed for a longitudinal 
field in , and a transverse field in Q . We focus in this 
paper on the synchrotron radiation by a moving quark 
in a longitudinal chromomagnetic field; this should be 
distinguished from the case studied in The author 
of [6( has considered the motion of fast fermion in an 
"electromagnetic" magnetic field not a chromomagnetic 
field. 

The synchrotron radiation of a photon from a relativis- 
tic electron was well studied in the last century 0, 
Generalizing the QED study to the QCD case has been 
done in [J] for a chromomagnetic field. It would have 
been the "end of a new chapter" but in [5] the path in- 
tegral approach was used to derive the synchrotron radi- 
ation of gluons by a parton in a transverse chromomag- 
netic field. The results found in Q seem to be radically 
different from Q; it was also argued that the QED case 
cannot be obtained from the QCD generalization of [|J. 
This motivates reconsidering the QCD generalization of 
the QED operator/source method used in [9j. We show 
that the QED case is simply obtained from our QCD gen- 
eralization. After clarifying the differences in the field 
configurations considered in [J] and [5|, it is possible to 
select the closest possible configuration to Hand make a 
comparison between our results and that of [5| . We show 



that very similar expressions are found. We analyze some 
of the possible inconsistencies that led to the erroneous 
results of 

The radiation in external magnetic fields is applica- 
ble in all Lorentz frames where H 2 — £ 2 > 0, where £ 
is the electric field. Hence we consider a weak field ap- 
proximation to avoid the "Klein catastrophe", that is, 
spontaneous pair creation by an electric field. Hence for 
\F llll F ,J ' v \ 1 / 2 <C m 2 /g, pair creation effects are negligible. 
Besides, for the semiclassical approximation to hold we 
have to have g 3 (F/m 2 )(E/m) <C 1; having large quark 
energy the only configuration to have a semiclassical ap- 
proximation is in the weak field regime, i.e., H <C m 2 /g. 
Therefore we consider the case where H <C m 2 / g. 

The paper is organized as follows. The coupling of 
quarks and gluons to the background field is first rep- 
resented in terms of color charges. The Schwinger et al 
method is then presented, the gluon synchrotron rate is 
calculated, and finally we compare the obtained results 
with different results in literature. 



II. CHARGE "REPRESENTATION" 

The main difference between the QED and QCD syn- 
chrotron is the interaction between the emitted gluon and 
the background chromomagnetic field. The non-Abelian 
generalization of the QED synchrotron radiation can be 
easily obtained if we characterize the coupling between 
gluons and the chromomagnetic background field through 
a fictive color charge. The Active charge can be explicitly 
obtained by starting from the QCD Lagrangian, decom- 
posing the gauge field into quantum fluctuations and a 
classical background field, then regrouping terms that 
couple to a given background color index. The coupling 
constant g is then found to be multiplied by a number 
which can then be interpreted as a color charge. The 
more formal approach is to use the Cartan subgroup 
of SU(3), then quarks and gluons can be characterized 
by two color charges. SU(3) is a unimodular group of 
3x3 Hermitian linearly independent matrices of deter- 
minant equal to 1. SU(3) has 3 2 — 1 = 8 generators A$ 
(i = 1 ■ • • 8). The maximum number of commuting gen- 
erators for SU(3) is 2. In the Gell-Mann fundamental 
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representation A3 and As commute hence T-L\ = A3/2 and 
% = A 8 /2 form the Cartan subgroup (they can be seen 
as the analogue of J3 for SU(2)) and are already diago- 
nal. Define the following combinations of the remaining 
six generators: 

E± a — Ai ± 1X2 E±p = Xq ± iX-j and E±~ ( — A4 ± iX$ 

these new matrices will play the role of raising and low- 
ering operators, they satisfy 

Hence the E% are eigenvectors of the ad(Hi) and ad(T~L2)- 
The weight £ of the adjoint representation are called the 
roots. They are the color charge of the gluon. A simple 
calculation, for the Gell-Mann representation, gives the 
gluon charges (£ (Hi), £(% 2 ) or (C^Qs)) 
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Quarks are found in the fundamental representation, so 
the eigenvalues of "Hi give the color charges of quarks. 
So quark "states" will be common eigenstates of ^1,^2 
with eigenvalues (93, q&)): 
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The color charge will always be multiplied by the strong 
interaction coupling constant g, so it is possible to absorb 
the "charge" in the coupling and define a new coupling 
9a — yQi f° r ghions and g l q — gq$. To simplify our 
calculation we will consider a chromomagnetic field with 
a given color index, either 3 or 8 but not combination of 
both color indices. 



III. THE MASS OPERATOR 

In the Schwinger et al method, the emission process 
can be seen through a modification of the Dirac equation 
of a spin 1/2 particle of mass m in the presence of an 
external field H . Besides the usual modification of the 
momentum operator (a tree level modification) to include 
the external field effect (i<9 M — > id — qA^ = 11^), the 
emission process can be seen as a modification of the 
mass by an operator M called the mass operator: 



M) * = 



So the mass operator is always acting on a physical spin 
1/2 real particle state. As we will see later, special atten- 
tion should be paid in some cases to its correct dispersion 
relation in the presence of an external field. 



The total decay rate T of a particle of mass m and 
energy E can be obtained from the mass operator M by 
an "optical theorem" -like relation 



T = -2— Im M 
E 



(2) 



where M is obtained from the one-loop diagram with full 
quark and gluon propagators in the presence of the field 
H . So the final expression will be valid to all orders in 
gH but not to all orders in a s . The decay rate in eq. @ 
is a number. M is understood as an average of the mass 
operator over different spin states: 



M 



uMu 



However we are usually interested in finding the prob- 
ability P of emitting a gluon with a given energy u. In 
such a case the radiation power for a given energy is re- 
lated to Im M by the integral equation 



(3) 



Hence to get P{oj) we have to write M in an integral 
form as we will see later. 



IV. EFFECTIVE PROPAGATORS 

Having a strong magnetic field, it is necessary to eval- 
uate the mass operator to all orders in the field H. We 
consider an "external" (classical) chromomagnetic field 
along the longitudinal direction (the third axis or the z 
axis): 



where the color index i is either 3 or 8. This field will lead 
to gauge field A that is linearly dependent on position. 

It is necessary to use effective quark and gluon propa- 
gators in the presence of H . 

The quark propagator in a field H can be written in a 
Fock-Schwinger proper time form as 



S(x, x') = <p q {x, x') 



d 4 p 

(27T)* 



exp[ip(x - x')]S(p) (4) 



where, the non-translation invariant, gauge dependent 
factor 4> is a Bohm-Aharonov-like phase 



4>q(x,x') = exp[igq 



(x)] 



(1) The "momentum" space quark propagator is 



s (j ) ) = i f -^-exp{-is 2 [ 
Jo cos z 
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1 ~E i 



(1-P)\\ 



tan z 9 



;(7-p)||> (5) 
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We have explicitly used the longitudinal and transverse 
decomposition appropriate for the choice of the field: (a • 

b)± = 0161+0262, ( a 'fr)|| = a>ob — The parameter z 
depends on the color charge of the considered quark z — 
gqhs2- The matrix S3 is the usual (z-)spin-projection. 
Graphically, in Feynman diagrams, the effective quark 
will have a blob to distinguish it from bare propagators. 

The nonzero coupling between the gluon and the back- 
ground leads to a modification of the bare propagator. 
The gluon effective propagator can be written in a form 
similar to the quark propagator 

G^(x,x ) = <f> Q (x,x') I —-^exp[ip(x-x)]G^ v (p) 



(6) 



where the gluon phase factor is again 

<p Q (x, x') = exp[igQ / {x)\ 

J X 

and the " momentum" representation gluon propagator is 



cosy 



dsi r , o tanw 
— exp{- lS i[-pf + — -pl]}E^ (7) 
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where y — gQhs\, and the tensor E^ v can be written 
(formally) in a compact form in terms of the tensor F^ v 
as 

Ep, = [exp(2gQFsi)}n„ (8) 
A more practical form is 

= fl|L + & cos(2y) - sin(2y) . (9) 

The metric tensors stand for the longitudinal and 
transverse spaces metric. The antisymmetric tensor A is 

A^ u = 5^\5 V 2 — dfisdui ■ 

It is easy to verify that in the limit of the vanishing mag- 
netic field E^v -> g^y . 



V. ONE-LOOP CONTRIBUTION 

The order a s contribution to the "mass operator" is 
given by the diagram shown in Fig. [TJ The coordinate 
space M is given then by 

j-^exp]ip(x-^)W(p) (10) 

where we have used the relation between the external 
quark charge q and the the charges of the internal quark 
q' and the gluon Q: q = q' + Q to write 

<j) q (x,x') = 4> q >{x,x')(j) Q (x,x') . 




p — k 

FIG. 1: The one-loop diagram contributing to the mass oper- 
ator in momentum space. Gluon and quark propagators are 
effective propagators in the presence of an external chromo- 
magnetic field. 



The "momentum" representation mass operator is 

M( P ) = -g 2 (T A ) ik (T A ) ki J -0^G,„(k)rS(p-k)r ■ 

(11) 

To get the un- integrated emission probability P(oj), we 
use the insertion used in [9(: 

00 00 

1= J dui J — exp{— i\(uj — k )} . 

—00 —00 

From now on our procedure will follow the steps of 
To avoid hindering the main result of this paper by tech- 
nical details we present the main steps and most of the 
derivation Appendix [X] 

The exact expression of the one-loop momentum rep- 
resentation M(p), after k and A integration, is 



00 00 



M(p) = -2-^--(T A ) ik {T A )u I du; I ds 



du 



5' cos y cos z 



9 rrr<A\ (rr A\ 

Hi 

exp(is<I>) \ me tz ^ 3 



(1 — u) tany 
yS' cos z 



+ [m - (p - u) lo + Pz (l - u)j 3 ] e^^^U) 

Few definitions are made: y = (1 — u)sgQh; z — usgq'h 

, .tan-y tan z 
6' = (1 - u) + u 

y z 

and the phase factor 

2 /-1 \ 2 / \ 1 z tan y 1 ,n 

<t> = —m +u{l — u)p» — u(l—u)— '- + {up — Loy . 

11 z y 

Note that M(p) at this stage is still in a matrix form. 



VI. FOURIER TRANSFORM: DISPERSION 
RELATION 

The expression of M(p) found in the previous section 
[Eq. (|12p] has to be Fourier transformed and then multi- 
plied by the gauge field depending phase <fi q (x,x'), with 
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the external quark charge q, to give M(x, x'). This can 
be done using the relations derived by Tsai in [lfj. For 
instance 



[x, x 



j4 

" P ip(x-x') iAavl iA 3 p% -iBp 2 , 



(2tt)< 



= cos Z(x\e lA ° n °e lA * n h lST1 ^ \x') (13) 

The n's are the generalized momentum operators in the 
presence of a magnetic field. This relation can be easily 
derived, as shown in Appendix [Bj if one notices that the 
right-hand side can be decomposed into three quantum 
non-relativistic propagators: a free particle of some fic- 
tive mass related to A Q living in the " o" one-dimensional 
space, a free particle living in the "3" one dimensional 
space, and finally a particle living in the transverse plane 
and under the action of an external magnetic field H per- 
pendicular to the plane. The three propagators are ex- 
actly known which give the above relation with Z = gqS 
and 

S = ta,n~ l (aqhB) 

gqh 

Similar relations are used for the product of exponential 
factors and p M to the left of Eq. (|T3"]l . The mass operator 
can then be obtained. This mass operator will be acting 
on physical quark states satisfying (11^7^ + m)ip = 0. 
Hence it is possible to do the following replacements 



IT 



nil + gqhY> 3 — m 2 



n 



p\\ 



where pn now represents the eigenvalues of ILi which is 
not affected by the longitudinal chromomagnetic field. 

Applying Tsai's relations and the external quark equa- 
tion of motion to Eq. ([12"]) the mass operator will be 



M 



cos Z 
du— exp 



5' cos y cos z 

„ (1 — u) tany 
+ cosZ- v 



me 



its / 

i(z-Z)S 3 



[(7-p)|i -m] 



yS' cos z 

+ [m - ( Po - uj) 1o + Pz (l- u) l3 ] e -*(*+«-a»)=» } (14) 



where Z = gqS, and 

S = — - tan" 
gqh 



tan z tan y 

su(l - u)— 

5' y 



defined as mentioned before from the Fourier transform. 
The new phase factor, which is a scalar now, is 

$' = u( Po - u) 2 - (1 - u)(up 3 - w 2 ) --(pl-pl- m 2 ) . 

(15) 



VII. GAMMA SUBSTITUTIONS 

The mass operator has a spinorial/matrix structure. 
It should be mentioned that we have used the exact dis- 
persion relation for the replacement of Oj_ in the phase 
factor, without solving the equation of motion. How- 
ever, the Dirac structure in front of the exponential fac- 
tor in (1141) depends on the gamma matrices in a non- 
trivial way. We use the same approximation done in 
the literature, which is the weak point of the method: 
the spinor structure of the incident quark is the same 
as that of the free particle, i.e., the Dirac spinor u s (p). 
So sandwiching the mass operator between u and u, the 
gamma matrices are replaced by scalar quantities. Sim- 
ple matrix calculation leads to the following replacement 

rules: 7 -> p /m, 73 p 3 /m, £3 (po - ^p^) ! 

2 

The -q — ± arc for the different helicities of the quark. 
The gamma replacement leads to a scalar mass operator 
that can be used to obtain P(u), 



M = -2 7 f T? (T A ) lfc (T A ) fcl / du I ds(-) I an 



cos Z 



5' cos y cos z 



(4tt) 2 



exp(zs$') {m(cos(z — Z) + cos(Z + z — 2y)) 



-V0 Po 



Po + m 



(sin(Z -z) + sin(Z + z- 2y)) 



cos Z 



(1 — u) tan y f p Q — p\ — m 
yS' cos z 



(p - oj)—cos(Z + z-2y) + ^(1 - u) cos(Z + z - 2y) 
m m 



+irj(p - u) 1 + 



Pi 



m(p + to) 



sin(Z + z - 2y) 



-ir/(l - «)^% (po ) sin(Z + z - 2y) 

Til \ p + TO ' 



Note that this expression has helicity information that 
can be explored in a way similar to [6|. 



A. The Abelian limit 

The results found in [9( can be easily obtained from 
our expression in Eq. (1161) . The non- Abelian character 
can be waived if we set the gluon charge Q to zero hence 
y — 0; remove the color factor T A T A in M; and consider 
the motion of the incident quark in the transverse plane, 
i.e., take p 3 = 0. Besides, we have to set the index of 
refraction n = 1 in [§]; hence the parameters 5, E, (3, 6' 
of [9( will be in our notation (for y = = p 3 ) 8 = 1, 
E^p , 5' -> 5', /3->- Z. 
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The above simplifications lead to 

,,2 



where £ = i — z, r = (z,p). For fixed £ the "transverse" 
wave function satisfies the equation 



L abenan 



(4tt) 2 



e?w / ds 



g?w exp(is<I>') 



{cos(Z — z) — sin(Z - 
m 

H — cos(Z + z) — in— sin(Z - 

mm m 



m 

z) 



dz 



pi 



2E 



z 



Besides, in the low mass limit or for m <C E the above 
equation simplifies to 



1 



AV2 



cos(Z 



17) 



whereas in @ A 1 / 2 = (5' cos 2)/ cos Z. Eq. (fTT)) coincides 
with equation (3.39 a) of 

B. Small s limit 

By analyzing Eq. it is clear that the s-integration 
is dominated by small a s region. The small s expansion 
of the phase factor gives 



$' « (up 



2 2 

m u 



- [gq'uh - gQh{\ - u)f & - [u 2 (l - u) 2 (p 2 - p 2 ) - m 2 ^18) 

The ^-integration can then be approximated by a Gaus- 
sian integration: 



duf(u)e 



is(up —u) 2 



/(«< 



— ) — 

Po ZlPo 



(19) 



In this approximation the variable u is replaced by the 
gluon energy fraction x = Lo/p - After the u- integration 
the phase factor will be 



dz 



pi 



2E 



(24) 



Which is a two-dimensional Schrodinger equation with a 
mass E. The solution to the above equation can be easily 
obtained, 



lp±-p 



-*J 



(25) 



which can be interpreted as a plane wave. The £ de- 
pendence emerges only via boundary conditions for the 
transverse wave function. So we have an evolution equa- 
tion along each line £ = const. 

If we consider now a parton of color charge Q in an 
external field H, represented by a gauge field A^, the 
equation of motion becomes 



dz 



(p - gQA) 2 



2E 



+ gQ(A°-A 3 )U. (26) 



The field configuration considered in Q is such that A = 
= A° and A3 = [H x ^3; hence the magnetic field is 
taken to be in the transverse plane. The wave function 
is then simply found to be 



s$' -> -m 2 sx 2 - — s 3 [x 2 (l - x) 2 (p 2 - pi) - m 2 } (20) 

where we have defined the vector 

f = gq'xB. - gQ(l - x)H 

which represents a fictive magnetic force on an interme- 
diate gluon-quark system. 

C. Zakharov method 

Before proceeding to an explicit comparison of our re- 
sult with the one found in [5[; we present Zakharov's 
method in its simplest "display." 

Consider first a fast scalar particle whose first quan- 
tized wave function satisfying the Klein-Gordon equation 



d 2 



d 2 



dt 2 dr 2 



= m 2 ip 



(21) 



For a fast particle moving along the z-axis with large 
energy E, the wave function can be written as 



ip(t, r) = exp {-iE£) ^{z, £, p) 



(27) 



But now the transverse momentum is a solution of the 
" classical equation of motion" 



dz 



This is applied to the incident quark as well as to the 
outgoing gluon and quark [Til ]. 

The final result of Q can be written as 



dP 
dLdx 



2^ 



dr \gi 



exp 



f 2 T 2 



<J2 



2fi 24/x 



(28) 



where p = Ex(l—x) (an effective mass), g\ = Ca s (l— x+ 
x 2 )/x (non-spin-flip vertex factor), gi = C'a s m 2 x 3 /2p 2 
(spin-flip vertex factor), e = m 2 x 2 + m 2 (l — x) 2 , and as 
in our approach the fictive force of a quark-gluon system 



(22) is f = iF,. - (1 - x)F g 
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Comparing a particle propagation in a longitudinal 
field to that of a transverse field configuration is not intu- 
itive. However, in an infinite (nonrealistic) medium, and 
for a static (non-propagating) magnetic field, and with a 
vanishing electric field, it is the propagation direction of 
the incident parton which gives the meaning of the lon- 
gitudinal/transverse directions. So if in our approach we 
take p z = we will be the closest to the configuration 
of [5| , where the transverse direction is mapped into the 
longitudinal direction in the following ways: 
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• Present configuration: a particle moving in the 
transverse plane under the action of a longitudinal 
field. 

• Zakharov's configuration: a fast particle moving 
along the longitudinal direction and a transverse 
magnetic field. 

It is closest but not the same, since in transverse 
momenta are not set to zero. 

So setting p z = 0, p = E and s — r/(2/i) in Eq. ([20)1 
we obtain the same phase given in Eq. (|28[) for a zero 
gluon mass. However, not surprisingly, the pref actor 
of our approximated expression does not match that of 
Eq. 0. 



D. Shuryak et al 

It is now clear that our agreement with Q leads to the 
conclusion that the results of [|| are " defective" . As the 
author of Q did not give enough details about their ap- 
proximation method it is hard to trace the exact source of 
inconsistency in Q ■ However we expect that their asym- 
metric treatment of the quark and gluon propagators 
leads to a nonsystematic approximation/expansion. For 
instance, the important Bohm-Aharonov phase, which 
plays an important role in our derivation, was absorbed 
in an unjustified way. Consequently the Active force f 
was badly approximated. 



VIII. CONCLUSIONS 

The generalization of QED synchrotron radiation to 
the QCD case was done with the minimal possible com- 
plications. It is shown that the QED case can be ob- 
tained as a special limit. The comparison with the path 
integral approach of Zakharov has shown consistent re- 
sults. This will hopefully "solve" the debate about gluon 
synchrotron radiation. 

It is possible now to extend our method to include 
nonzero gluon mass (thermal mass) and combinations of 
longitudinal and transverse chromoelectric and chromo- 
magnetic fields with an arbitrary color index. This is 
under current investigation. 



Appendix A: Schwinger et al procedure 

Our derivation follows the recipe of @. We present in 
what follows the main steps needed to get equation (|12p . 

The first step to follow is the change of variables: si — 
(1 — u)s, $2 = us, where < u < 1, and s is positive. 
The A and k integration is based on Gaussian integration 
and saddle point approximation. Hence the exponential 
factor is factored out in the expression of M; this can be 
written as 



M(p) = J du 
o 



dXdoj J dsM D 
o 



</'/,- / d\dx I tisM, <«xp j is(x ) 

(Al) 

where M Dirac contains the Dirac structure of the self- 
energy diagram and all the non-exponential terms. The 
function \ can be split into a momentum dependent part 
Xi and a A dependent function \ a , \ = Xo + Xi, such 
that 



Xi = [k - (up 
while 

Xo 



2s 



-/ / , " tan z _ 
-(k 3 -up 3 ) -6 [kj_ J^—P^ 



4 



up Xs 



where 5' and $ are as defined before. 

We perform now the k integration using the Gaussian 
form(s) 



d 4 k 
i 



exp(isxi) 



1, k a , k 3 ,k± 



A u tan z _ 

1, upo - 7T,up 3 , — p± 

2s o'z 



(A2) 



S>(4ir) 2 s 2 
For the A integration 

/ 2~tt exp ^ SXo ~ A l 

( — ) V2 e »(*+(«P--« ) 2 ) [i, 2s(u Po - w)l (A3) 



The remaining, lengthy but straightforward, step needed 
to get equation (TT2"j) is the Dirac structure simplifica- 
tion / contractions. 
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Appendix B: Tsai transformation to coordinate 
space 



straight line trajectory. Hence the four-dimensional gen- 
eralization will be 



We give a simple derivation for equation (113[) which is 
different from the initial method of Tsai [1C§ . 

Consider a non-relativistic particle of mass m and 
charge g in a uniform magnetic held H . The general ex- 
pression of the non-relativistic propagator in coordinate 
space is 



K(r, r, f a , 0) = (f| exp {-iH^\ \r Q ) 



G 



m \ 3 / 2 gHT/2m 



tilths J sm(gHT)/2m 
where % is the particle Hamiltonian 

t2 m 



exp 



iS r 



(Bl) 



H 



51. _ 

2m 2m 



with II the generalized momentum. The classical action 
S c i in the presence of a magnetic field is 

5*d = ^-{~ [(x - x Q ) 2 + (y - y ) 2 ] cot( 5 i/r/2m) 

+ (x y - y D x)} + ^-(z- z ) 2 (B2) 

ZT 

The non-translational invariant term (x y — y a x) in the 
action gives the Bohm-Aharonov phase, if we assume a 



K(x tl ,T,x' tl ,0) = $(x,x') 



exp ( igH cot ) exp | -i - " "' ) (B3) 



\ 47tt y sin z 
.(^||-^|) 2> 



4t 



where z = g-ffr. It is now possible to write the propaga- 
tor K in momentum space. We use the Fourier transform 



' ;t / . (x_l — a;', ) 2 
■ exp igii cot z — — 



it sin z 



f d 2 p±_ . / . tanz 2 V-rws 

/ exp (-ip_L • (xj. - Xj_)) exp -it p ± (B4) 

J cosz \ z y 

This is similarly done for the longitudinal part (||) of the 
propagator. Hence 



K(x^,t, a;'", 0) = (x| exp (^n 2 ) |x') 



$(x,x') 



j,P- [x— x ) 



,»TP|| 



cosz 



where we have set m = 1/2 and ft = 1 in the nonrelativis- 
tic formula. This Formula is the same as equation ([13")) 
with appropriate relabeling of field and time parameters. 
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